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ABSTRACT 

We construct catalogues of present superstructures that, according to a hCDM scenario, will 
evolve into isolated, virialized structures in the future. We use a smoothed luminosity density 
map derived from galaxies in SDSS-DR7 data and separate high luminosity density peaks. 
The luminosity density map is obtained from a volume-limited sample of galaxies in the 
spectroscopic galaxy catalogue, within the SDSS-DR7 footprint area and in the redshift range 
0.04 < z < 0.12. Other two samples are constructed for calibration and testing purposes, up to 
z = 0.10 and z - 0.15. The luminosity of each galaxy is spread using an Epanechnikov kernel 
of 8Mpc/h radius, and the map is constructed on a 1 Mpc/h cubic cells grid. Future virialized 
structures (FVS) are identified as regions with overdensity above a g iven threshold, calib rated 
using a AC DM numerical simulation, and the criteria presented by lDiinner et alJ (l2006l) . We 
assume a constant mass-to-luminosity ratio and impose the further condition of a minimum 
luminosity of lO'^L©. According to our calibrations with a numerical simulation, these criteria 
lead to a negligible contamination by less overdense (non FVS) superstructures. We present a 
catalogue of superstructures in the SDSS-DR7 area within redshift 0.04 < z < 0.12 and test 
the reliability of our method by studying different subsamples as well as a mock catalogue. We 
compute the luminosity and volume distributions of the superstructures finding that about 10% 
of the luminosity (mass) will end up in future virialized structures. The fraction of groups and 
X-ray clusters in these superstructures is higher for groups/clusters of higher mass, suggesting 
that future cluster mergers will involve the most massive systems. We also analyse known 
structures in the present Universe and compare with our catalogue of FVS. 

Key words: large scale structure of the universe - statistics - data analysis 



1 INTRODUCTION 

Structures in the universe are the result of a hierarchical pro- 
cess of accretion dominated, in almost all scales and during al- 
most their entire history, by gravity. The first attempt to investi- 
gate the evidence for the existenc e of the so call ed second-order 
clusters of galaxies was made by IShaplevI l fT96lh . although there 
had been previou s suggestions of the presence of large structures. 
IShaplev & AmesI (IT930) suggested the presence of a large galaxy 
system in the Coma- Virgo region. Also, conglomerates composed 
by several clusters of galaxies observed in plates of th e Lick As- 
trographic Su rvey were reported by Shane & Wirtane^ ( 1 19541) and 
IShanef ( ll956h . By analyzing the distribution of rich clusters identi- 
fied on the National Geographic Society-Palomar Observatory Sky 
Survey, Abell pointed out that the clusters did not appear to be dis- 
tributed randomly over the sky, but forming associations of matter 
on greater scales. Different authors adopted a variety of methods 
to search for groups of clusters and the term supercluster became 
widely used. Abell used counts in cells to statistically test the distri- 
bution of clusters, and endorsed the existence of superclusters of or- 
der 50 hr^ Mpc in size. The first attempts to study superclusters on 



a statistical basis were performed by lin king Abell cluster positions 
jZucca et al ]| 19931 : lEinastoetalJI 19971) . Later, the accomplishment 
of wide-area surveys of galaxies with s pectroscopic follow u 
such as Las Campanas Redshift Survey ('Shectman et al."l99i 
the 2- degree Field Galaxy Redshift Survey (2dFGRS, Colless et al 



2001') and the Sloan Digital Sky Survey (SDSS. IStoughton et al 



2003), allowed for the identificatio n of superclus t ers dir ectly from 



the large-scale galaxy distribution. lEinasto et alj ( l2007h identified 
superclust ers in the 2d FGRS using a density field method, and re- 
cently Costa:;Duarte_elLal. (2010), studied the morphology of su- 
perclusters of galaxies in the SDSS. The largest catalogue of su- 
perclusters has been constructed by iLiivamagi et al.l ( |2010|) . who 
implemented the density field method on the SDSS-DR7 main and 
LRG samples. In all cases, superclusters are operationally defined 
as objects within a region of positive galaxy density contrast and 
thus are subject to a certain degree of arbitrariness in the parameter 
selection. 

Within the AC DM Concordance Cosmological Model, an ac- 
celerated expansion dominates the present and future dynamics of 
the universe and thus determines the nature of gravitationally bound 
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structures. Therefore, an alternative definition of these large-scale 
structures may be derived from the properties of overdense regions 
in the present-day universe, that will become bound and virialized 
structures in the future. Thus, under the assumption that luminosity 
is a somewhat unbiased tracer of mass on large scales, the inte- 
grated luminosity density of galaxies is commonly used as an indi- 
cator of mass density. 

The cosmological evolution of the large scale structure has 
implications on the spatial distribution, frequency and properties of 
superclusters and the galaxies they contain. Due to this relationship, 
superclusters can be used as cosmological probes, and their study 
is oriented to constrain models and describe the formation of su- 
perstructures on a cosmological scale. Supercluster properties have 
been used to discriminate between cosmological models, favoring 
the s t andard cosmologic al model in most cases jBasilakos et al.l 
l200lHPeacock et al.1 '2001). On the other hand, there are claims of 
structures that a re too massive or fo rmed too early according to the 
standard model teaugh et alj2004h . 

First attempts to determine the characteristic scales of spatial 
inhom ogene ities in the un i verse were made by Br oadhurst et al.l 
and lEinasto et al] ( Il994) . By studying the distribution 
of rich clus t ers of galaxies from the Abell-ACO catalogue, 
lEinasto et j 19941) confirmed a 110 - lAQlr^Mpc scale in the 
supercluster-void netwo rk. Given the obse rved regular pattern of 
superclusters and voids, iFrisch et al.l j 19951) investigated the prop- 
erties of the initial power spectrum giving rise to these large-scale 
fluctuations. The authors found that the supercluster-supervoid net- 
work forms in a very early stage of the evolution of the Universe 
from large-scale density fluctuations, and are defined by the scale 
of the maximum of the power spectrum. 

The large scale structure of the universe is often described 
as a supercluster-void network JShandarin et al. 2004), and su- 
perclusters are closely re lated to fila ments i Gonzalez & Padillal 
20101; iMurphv et al.l I20T0I) and void s jEinasto et alj|l986l ll997F 



Park & Lej|2007l: IPlaten et alj|2008l; Ilckelll984l) . Such filaments 



may als o play an importan t role in the process of structure for- 
mation. IPorter et alj ( |2008|) find that star formation is signifi- 
cantly enhanced when galaxies fall into clusters along superclus- 
ter filaments. S uperclusters have a l so been studied as hosts of 
Lytt absorbers ( IStockeetalJll995l : | Penton et al. I l2003) and are 
known to p roduce signatures in the Cosmic Backg round Radiation 
jGranett e t al. 2009; Dolag et al. 2005; Genova-Sant os et alj|2008l ; 
iFlores-Cacho et al.„2009 ; .Genova-Santos et al..2010l) . 

The aim of this work is to present and analyse catalogues of 
superstructures that will evolve into v i rialize d systems. Using the 
theoretical framework o f iDunneretalltoel) in a AC DM scenario 
and by calibration with numerical simulations we analyse volume 
limited samples of galaxies from the SDSS-DR7. 

This paper is organized as follows. In section Section |2] we 
describe the data used in the computation of the density field. The 
methodological description of the procedures implemented to ob- 
tain the catalogue of future virialized structures (hereafter FVS) is 
addressed in Section [3] along with a brief description of the pre- 
vious methodology employed in the search of superclusters. The 
catalogue of FVS is presented in Section|4] where we also discuss 
some of their properties. In Section |5] we compare known super- 
clusters with our identified structures, and in section Section|6]we 
study the frequency of clusters and groups of galaxies in future 
virialized structures. Our conclusions are summarized in Section 
[8] Throughout this paper, we adopt a concordance cosmological 
model (Da = 0.75, Sl„,a,ter = 0.25) in the calculation of distances. 



2 DATA AND SAMPLES 

Since FVS are extended regions of several tens of Mpc in size, a 
large volume of space has to be surveyed in order to sample a rep- 
resentative population of these objects. The Sloan Digital Sky Sur- 
vey (SDSS, Stoughton et al. 2002) is the largest photometric and 
spectroscopic survey carried out so far, covering an area of almost 
a quarter of the sky with a limiting magnitude that makes the con- 
struction of complete samples of several hundred Mpcs possible. 
Several catalogues have been made public by successive releases 
since Early Data Release ( Stoughton et al. 2002 ). The latest release 
of the spectroscopic catalogue (DRV. lAbazaiian et al. I I2OO9I com- 
prises 929 555 galaxy spectra within a footprint area of 9380 sq. 
deg. The limiting apparent magnitude for the s pectroscopic cata- 
logue in the r-band is 17.77 dStrauss et ai] |2002l). although we use 
a more conservative limit of 17.5 to ensure completeness. We also 
limit the sample to galaxies fainter than r = 14.5, since saturation 
effects in the photometric pipeline does not secure completeness 
below that limit. These limits were adopted taking into account the 
analysis of image quality and efficiency detection of the SDSsQ 

Given the dependence of the luminosity density field on the 
sample of galaxies, we define three complete samples, with dif- 
ferent cuts in luminosity (Table [TJ. The closest sample comprises 
galaxies brighter than M,- = -20.5 up to z = 0.1. This is the sam- 
ple with the faintest luminosity limit, and it will be mainly used 
to explore the effects of the luminosity cut in the detection of 
superstructures. This sample will be referenced as SI, and their 
characteristics are summarized in Table [T] A larger volume lim- 
ited sample, referenced as S3, contains all galaxies brighter than 
M,- = -21.0 to z < 0.15. An intermediate redshift sample in the 
range 0.04 < z < 0.12 comprises 89513 galaxies. This sample, S2, 
containing galaxies with M,- < -20.47 will be analised into detail 
in Sections|4l|5]and[6l 

In order to test procedur es and results, we use the Mil- 
lennium numerical simulation ( Springel et al.ll2003) of a AC DM 
concordance model, performed on a cubic box of 500 Mpc 
side. A semi-ana lytical model of galaxy formation (GALFORM, 
iBower et al. II2OO8) collects information from the merger trees ex- 
tracted from the simulation, and generates a population of galaxies 
within the simulation box. We use the semi-analytic galaxy cata- 
logue in the full box of the simulation to test and set the parameters 
involved in the identification of FVS. We constructed a mock cat- 
alogue from the semi-analytic galaxies following the geometry of 
the SDSS-DR7 footprint area and reproducing the dilution in the 
number of galaxies with redshift. In order to test the effect of pe- 
culiar velocities in the identification of FVS, we defined a sample 
of galaxies in real-space Mr,,, and a sample of galaxies in redshift 
space Mzsp, also described in Table[T] 

We use AB magnitudes and apply k-He correctio ns in the rest 
frame at z = 0.1, as defined in telanton et alj|20()3l) . From these 
magnitudes, we calculate the luminosity of each galaxy, in the r- 
band, as L = Lq x io.<o-+x(*'o-m,)) ^^^^^^ luminosity density 
Plum of each sample is computed with this luminosity, and is also 
listed in Table [T] We will also study the properties of groups and 
clusters of galaxies within FVS. To this end, we use a catalogue o f 
galaxy groups in the SDSS-DR7 compiled by'Zapata et al.' ("2009"). 
This catalogue is obtained using an adaptative Friends-of-Friends 
method. Virial masses are computed from the line of sight velocity 
dispersion of galaxies in eac h group (cr,,) and the virial radius (/?,.,>), 
using the virial theorem (see lZapata et al.ll2009h : 



http://www.sdss.org/dr7/products/general/target_quality.html 
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Sample 




D,„„,r/i"' Mpcl 

lllCl.\ L" irj ^^L. J 


M^'"' 


Volume [10'(/i"' Mpcf] 




p,,,„,\W^Lr^|Mn^'^ 


F 


CoiTected pi„,„[10^ Lq/M pc^] 


SI 


0.10 


293.92 


-20.05 


1.85 


94271 


0.80 


2.11 


1.68 


S2 


0.12 


351.34 


-20.47 


3.17 


89513 


0.58 


2.98 


1.73 


S3 


0.15 


436.55 


-21.00 


6.01 


62344 


0.29 


5.66 


1.64 


S2c 


0.10 


293.92 


-20.47 


1.85 


51188 


0.56 


2.98 


1.73 


S3c 


0.10 


293.92 


-21.00 


1.85 


17507 


0.27 


5.66 


1.64 


Mrs,, 


0.12 


351.34 


-20.47 


3.17 


106604 


0.75 


2.98 


2.23 




0.12 


351.34 


-20.47 


3.17 


106722 


0.75 


2.98 


2.23 



Table 1. Galaxy samples in the SDSS-DR7. In all cases: z > 0.04, and the apparent magnitude in the r-band is in the range 14.5 ^ r ^ 17.5. The mean 
luminosity density p;„„, is computed using volume limited samples, each containing Ngai galaxies. The connection factor F (Eq. |6) and the resulting mean 
luminosity p/„m of each sample are indicated in the table. 
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(1) 



where i?,-i> is estimated as in iMerchan & Zandivare j j2005t) . The 
catalogue comprises 83784 groups with at least 4 members, and is 
limited to redshift z < 0. 12. 

We have also a nalysed a catalogue of X-ray selected clusters 
jPopesso et alj|2004l) in order to study the presence of large gravi- 
tational potential wells in FVS. This sample comprises 114 clusters 
with X-ray emission, and is based on the ROSAT All Sky Survey 
(RASS) and the Sloan Digital Sky Survey (SDSS). The total lu- 
minosity in the ROSAT band is avail able for each clust er, and is 
used as a proxy for the cluster mass tevkoffet alj|2008l) . We use 
this luminosity in Section [6] to explore the variations with X-ray 
luminosity of the fraction of clusters belonging to FVS. 



3 METHOD 

3.1 Previous analyses 

Several procedures have been proposed to construct a reliable cata- 
log of superstructures, although they are based on different criteria 
and methods. The main issue related to the comparison of super- 
structure catalogues is that the definition of superstruc tures depends 
on the specific choice of a number of free parameters. lEinasto et al.l 
use a percolation algorithm to link clusters of galaxies from 
the Abell-ACO catalogue, and derive a catalogue of 220 superclus- 
ters of rich clusters within z = 0.12. The authors percolate clusters 
lying within a given radius from each other, defining these systems 
of clusters as superclusters. They find that a percolation radius of 
24 h^'Mpc is convenient to detect the largest, but still relatively 
isolated, systems of clusters. The main results from the superclus- 
ter catalogue are not, according to th e authors, very sensit ive to the 
exact value of the percolation radius. Isinasto et al] i200lh updated 
the catalogue with a newer version of the Abell catalogue, incorpo- 
rating to the analysis a sample of X-ray clusters; this also allowed 
them to compare superclusters derived from different samples of 
clusters. They find that both types of clusters generate superstruc- 
tures that represent the large scale structure in a similar way. Also, 
they find that X-ray clusters not belonging to superclusters sur- 
round the Southern and Northern Local supervoid, or are located in 
filaments between superclusters. Although the authors find a strong 
signal indicating that the fraction of X-ray clusters in superclusters 
increases with supercluster richness, there is no correlation between 
the X-ray luminosity of clusters and their host sup ercluster rich- 
ness, quantified as the number of member clusters. lEinasto et al.l 
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Figure 1. A slice of the luminosity density map for the S2 sample. The 
cartesian system used is associated with the equatorial coordinate system, 
with the z-axis towards the north celestial pole. The smoothing kernel has 
an Epanechnikov shape, with a size of 8 ft" 'Mpc 



( l2007h return to the problem of defining the largest i solated struc- 
tures in the universe, and use data from the 2dFGRS jCoUess et al.l 
2001) to assemble a catalogue of superclusters. The authors find 
that the most effective method to perform a supercluster search is 
the density field method which consists in obtaining a smoothed 
luminosity density field in redshift-space from the galaxy cata- 
logue. The luminosity, on large scales, is supposed to follow the 
distribution of matter, provided that a conv enient smoothing kernel 
is used. In their work, lEinasto et alj | |2007!) use an Epanechnikov 
kernel, and assert that the best results in the identifi cation method 
are ob tained when a kernel size of 8 /i"' Mpc is used. lEinasto et al] 
( l2007t) set the threshold parame ter by maximi zing the number of 
large superclusters. Costa-Duart e et al] ( 1201 Ot) apply the density 
field method to the SDSS-DR7 and use two samples of structures, 
selected using different overdensity thresholds. The authors claim 
that there is no natural value for the threshold density, and define a 
sample that maximizes the number of structures, and a sample with 
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the parameters tuned so that the largest superclusters present an ex- 
tension of a 120 h"' Mpc. The authors perform an analysis of the 
shapes of superclusters using Minkowski functionals, finding that 
filamentary structures tend to be richer, larger and more luminous 
than pancakes. They also use a semi-analytic catalogue of galaxies 
derived from the Millennium simulation to test the method and con- 
clude that the morphological cla ssification is not biased by peculiar 
velocities. Liivama gi et alj j2010 ) use a similar approach, but per- 
form the smoothing with a B3-spline kernel of radius of 8/i' Mpc, 
obtaining catalogues of superstructures with similar properties in 
the SDSS-DR7 and in the Millennium simulation. 



3.2 Present approach 

Given that not all large-scale structures are virialized at the present 
time, the setting of identification parameters are subject to certain 
degree of arbitrariness. A physically motivated thres hold on mass 
overde nsity was explored in AC DM simulations bv iDunner et al.l 
( l2006l) . According to these authors, it is possible to define a cri- 
terion to isolate, using three-dimensional data, overdensity regions 
enclosed by a spherical shell that will evolve into virialized sys- 
tems. By the application of the spherical collapse model, the mean 
mass density enclosed by the last bound shell of a structure must 
satisfy: 

rrmass 

^ =7.88, 

f^bck 



(2) 



where p")°;/ is the mean mass density enclosed by the critical shell 
(the shell that maximizes the potential energy), and p)"^°:'' is the 
mean density of the background. In observational catalogues there 
is not an accurate estimation of the mass density field. However, 
given that at large scales the mass-luminosity ratio is nearly con- 
stant, we can apply a similar criterion to the luminosity map to 
derive structures with an appropriate mass overdensity. 

In the next subsections we describe the identification method 
implemented on the SDSS-DR7 spectroscopic data. We first define 
the volume covered by the sample using a three-dimensional mask 
and construct a luminosity density map with a 1 (h^' Mpcf cell 
resolution. We then use a percolation method based on the search of 
high density peaks on the three-dimensional smoothed map. These 
overdensities are the basis of the superstructure catalogue. 



3.2.1 Luminosity Density Map 

We construct smooth luminosity density maps from our complete 
samples of galaxies, in the region defined by a three-dimensional 
mask. This mask, which represents an approximation of the geom- 
etry of the SDSS-DR7 galaxy catalogue in both angular coordi- 
nates and redshift, is built using cubic cells of 1 Mpc side. We 
start from a pixelized representation of the central region within 
the catalogue footprint area. To this en d we use an angula r mask 
obtained using the software HEALPIX dOorski et al.ll2005 f) with a 
resolution parameter A'„rf(, = 512 that splits the whole sphere in 
3 145 728 equal area pixels. A cubic cell is part of the 3D mask if 
at least a fraction of its area includes part of the solid angle sub- 
tended by the angular mask, and its radial position lies within the 
redshift range defined for a given sample. We define the volume of 
a cell within the three-dimensional mask as V'^^n = V x fceii, where 
V = (l/i"' Mpc)^ and /„// measures the fraction of the volume V 
into the mask. We compute the fraction /„/; by implementing a 
Monte Carlo procedure. Accordingly, cells that are completely con- 
tained within the geometry of the catalogue have f^eii = 1, while 
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Figure 2. Distribution of the mass overdensity of the superstructures identi- 
fied in the simulation box, for different values of the luminosity overdensity 
parameter (Dt = 4.5, 5.5 and 9, in the top left, top right and bottom left pan- 
els respectively). The dashed vertical lines indicate the critic al mass over- 
density that a region must have to be virialized, according to lPiinner et alJ 
Ham. Dotted histograms show the distributions for the total sample of 
superstructures, and black continuous histograms show the con'esponding 
distributions for structures with total luminosity greater than IO'^Lq. Bot- 
tom right: comparison between the distributions of mass overdensity for the 
three values of Dr. 



external cells have = 0. The factor f^^n allow us to correct for 
border effects, in our subsequent analysis we will use only cells 
with fcfii > 0.5 i.e., only cells with at least 50% of their volume in- 
side the mask. The continuous luminosity density map can be con- 
structed by smoothing the galaxy distribution within the 3D mask. 
A standard procedure to accomplish this consists in using a kernel 
function to convolve the discrete positions of galaxies and spread 
their luminosity. The resulting density field is then represented at a 
resolution given by the cell size, and corrected by the weight fac- 
tor /„;;. The result of the smoothing depends on the shape and size 
of the chosen kernel function. Following previous analysis in the 
literature (Ei nasto et a l. 2007 ; Costa-Duarte et al . 2010). we use an 
Epanechnikov kernel of ro = 8/!" 'Mpc size, which gives the con- 
tribution at position r from a source that is located at R: 



k(r - R) = — 

4ro 



|r-R| 

ro 



(3) 



An Epanechnikov kernel is more suitable for this analysis since its 
shape resembles that of a Gaussian, but it avoids excessive smooth- 
ing. The luminosity density estimate within a cell is then: 



Peel I - LcelllVcei 



(4) 



where V^^,, = V x f^eii is the volume of the cell and L„;; is the sum 
of the contributions to the luminosity L^;, from nearby galaxies: 



^ ^ cell 



k(r - Ri)dr. 



(5) 
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Figure 3. Completeness and contamination of the FVS catalogue as a func- 
tion of Dj. The chosen value of Dj = 5.5 for the compilation of the cata- 
logue of FVS shows the best balance between low contamination and high 
completeness. 



A slice of the map is shown in Fig.[T] where the large scale structure 
of the luminosity distribution can be appreciated. Concentrations of 
galaxies such as rich groups or clusters are characterized by lumi- 
nosity density peaks, and are surrounded by low density regions. 



3.3 Identification of structures 

We search for large isolated regions above a given density thresh- 
old; these are candidates to FVS. The luminosity density, however, 
depends on the limiting magnitude characterizing a given volume 
limited sample, which determines the number and luminosity of 
galaxies contributing to the overall density estimate. For example, 
deeper volume limited samples include only brighter galaxies, thus 
the luminosities of the cells will be systematically lower as the up- 
per redshift limit of the sample increases. We use the luminosity 
overdensity instead of the total luminosity density value to charac- 
terize the peaks in the luminosity distribution, but still, since we 
study three samples with different limiting redshifts, a correction 
factor must be implemented in order to make the results as compa- 
rable as possible. This factor can be obtained by assuming a uni- 
versal luminosity function for galaxies <I>(L), which allow us to ac- 
count for the luminosity below the limiting value characterizing 
a given sample. Since the underestimation of total luminosity de- 
pends only on the absolute magnitude cut M;,„ of the sample, an 
homogeneous correction factor F for the entire luminosity density 
map can be defined as: 



X°° m{L)dL 



(6) 



We ad opted the luminosity function presented by iBlanton et al.l 
( l2003l) . This correction ranges from 2 to < 6 for samples with 
limiting absolute magnitudes down to = -21. In Table [T] the 



correction factor F and the corrected mean luminosity density are 
listed for each sample. 

We define superstructures by linking overdense cells using a 
Friends-of-Friends algorithm that connects overdense cells sharing 
at least one common vertex or side. To this end, we use a thresh- 
old luminosity overdensity criterion pium-ceii > Djpu,,,,. Since the 
properties of the catalogue of superstructures may be affected by 
the adopted value of the luminosity overdensity parameter Dt, it is 
fundamental to study this issue into further detail. 



3.4 Calibrating the method with numerical simulations 

We use numerical simulations to test the ability of a luminosity- 
based algorithm to derive structures that will evolve into virialized 
systems in the future. We have also used the simulations to analyse 
whether superstructures identified in redshift-space with total lu- 
minosities above L = 10'^ Lq, correspond to systems in real-space 
above this threshold. This is an important issue to be checked since 
the observations provide redshift-space data whereas the FVS cri- 
terion requires real-space information. In Fig.|2]we show the dis- 
tributions of mass overdensity (only known in the simulations) for 
the sample of structures that result from the identification algorithm 
for three different values of Dt, and the corresponding distributions 
for superstructures with a luminosity above the threshold 10.'^ Lq. 
This adopted luminosity threshold is approximately the maximum 
observed value for a rich cluster of galaxies; by imposing this lu- 
minosity cut-off we exclude individual clusters from our FVS cata- 
logue. The bottom right panel of Fig.|2]shows the mass overdensity 
of superstructures with the luminosity cut-off for three values in 
luminosity overdensity parameter. The vertical dashed line in this 
figure shows the minimum mass over density for a superstr ucture to 
be virialized in the future according to Dtin ner et alj 00061) . As can 
be seen, the contamination in the catalogue depends on the choice 
of the luminosity overdensity parameter. In order to select the most 
convenient value of this parameter, we explored the contamination 
and completeness of superstructures for different values of Dj. We 
define the completeness as the ratio of the number of superstruc- 
tures resulting from a given value of Dj and maximum number of 
superstructures obtained within the explored Dj range. Similarly, 
the contamination is defined as the fraction of identified superstruc- 
tures that do not satisfy the future virialized criteria. In Fig. [3] we 
show the completeness and contamination parameters for the range 
4 < Dt- < 9. By inspection of this figure, we chose a luminos- 
ity overdensity threshold Dj = 5.5. This adopted value provides a 
suitable compromise of high completeness and low contamination. 



3.5 Redshift-space vs. real-space analysis 

In this subsection we analyse the effects of redshift-space distor- 
tions induced by peculiar velocities on the identification of FVS, 
using the mock catalogue. We considered the samples and 
Mzsp (Table [Xll which use the real-space and redshift-space po- 
sition of mock galaxies, respectively, and are defined using the 
same redshift limits and limiting absolute magnitude than sample 
S2. This allows to study FVS identified in real-space and compare 
them with those selected in redshift-space. We have computed the 
distributions of luminosity and volume derived for the FVS in these 
mock samples. A good agreement can be appreciated between the 
luminosity distributions (Fig. |4(a)^ , indicating that the identifica- 
tion procedure delivers reliable estimations of the total luminosity. 
It can also be inferred from the volume distributions, shown in Fig. 
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Figure 4. Histograms representing the luminosity and volume distributions of Future Virialized Structures. Solid lines correspond to the sample of FVS 
identified in sample Mzsp (redshift-space mock, see Section|2), and dotted lines coirespond to FVS identified in sample Mr,,, (real-space mock). 



|4(b)[ that the volume of FVS are not strongly affected by peculiar 
velocities. 

We also analysed the correspondence between structures in 
real and redshift-space in order to estimate the contamination and 
completeness of the FVS catalogue. The first point we notice is 
the lower number of FVS identified in the redshift-space sample in 
comparison to the real-space sample, in the same volume. About 
20 per cent of FVS identified in real-space are lost when analyzing 
redshift-space data. This is the major drawback of redshift-space 
FVS identification. On the other hand, only 2 per cent of FVS iden- 
tified in redshift-space are not real-space FVS. This indicates that 
the identification of FVS in observational catalogues is not likely 
to include fake FVS structures. We also found that approximately 3 
per cent of FVS in the mock catalogue are associated to more than a 
single FVS in real-space. A similar percentage of real-space FVS 
are identified as multiple FVS in redshift-space. 



4 CATALOGUE OF FVS 

Taking into account the results of the previous sections,we have 
adopted a luminosity overdensity threshold of Dj = 5.5 and in 
order to avoid the inclusion of spurious systems, a lower luminosity 
limit > IO'^Lq- In Table[2]we summarize some characteristics 
of the identified FVS. 

Once the superstructures have been identified, it is important 
to assess their fundamental properties such as total volume, number 
of galaxies above a given luminosity, the total luminosity, shape 
parameters, etc. Given the small volume of sample SI, and the large 
luminosity correction factor in sample S3, we provide an analysis 
of sample S2 in what follows. 



Sample 


Nfvs 


Prol 


Finn, 


glxsi„fvs 


SI 


67 


1.08% 


10.85% 


9707 


S2 


150 


1.26% 


13.54% 


11394 


S3 


412 


1.66% 


20.61% 


11682 


Mr.j, 


227 


1.62% 


18.87% 


19265 


Mz,,, 


181 


1.35% 


15.14% 


15368 



Table 2. Main results obtained for the samples of identified FVS. For each 
sample, we show the number of future viriaHzed structures Nfvs • the per- 
centage of volume occupied by FVS F,.o/. the percentage of luminosity of 
galaxies within FVS -F;„„, and the total number of galaxies within FVS 
g!xsi„fvs ■ 



4.1 Analysis of the FVS catalogue (sample S2) 



In Fig. |5(a)| we show the distribution of luminosities of the super- 
structures we obtain for the SDSS-DR7. FVS luminosities vary be- 
twe en IO'^Lq and ~ 10^* Lg, in agreement with Einasto et al. ( 20061) 
and lCosta-Duarte et al. However. lEinasto et al.i (2006.) find 

a lack of luminous superclusters in numerical simulations com- 
pared to superclusters identified in observational catalogues. In a 
forthcoming paper (Luparello et al., in preparation) we will anal- 
yse this issue in more detail. 

Since the superstructures are obtained from a discrete density 
map, the volume can be directly calculated as the sum of the vol- 
umes of all the cells that belong to a given system. As all cells have 
the same volume, equal to 1 {hr^ Mpcf , the number of cells that 
form a structure is directly proportional to the volume. As it can be 
appreciated in Fig. |5(b)| the volume of the superstructures ranges 
between 10^(/!"'Mpc)' and \^{hr^ Mpc)^ . We also compute the 
fraction of the total volume that is in FVS, finding that FVS repre- 
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Figure 5. Luminosity and volume distributions of structures for samples SI (dashed line), S2 (solid line) and S3 (dotted line). The insets show the average 
voluine and luininosity as a function of the redshift, in redshift bins with equal numbers of superstructures, with filled circles, empty circles and cruxes, 
respectively for samples SI, S2 and S3. EiTor bars on the histograms indicate Poisson uncertainty. In the inset box, the uncertainty bars correspond to the 25 
per cent and 75 per cent percentiles. 



sents only L26 per cent of the volume covered by the sample S2. 
lEinasto et alj ( |2007|) found that the volume covered by superclus- 
ters in the 2dFGRS represents 3.2 per cent of the Northern region 
of the catalogue and 3.5 per cent in the Southern region. Our lower 
percentages can be explained by the fact that we use a higher value 
of Dj and impose a limit of at least lO'^Lg in the luminosity of the 
superstructures. 

Previous analyses of numerical simulations have shown that 
the shapes of halos are consistent with an increasing elongation 
as the mass inc reases. This is also detected in observational data 
jPaz et alj20"0^ . where the most massive structures tend to be pro- 
late. We notice that the shapes of FVS are somewhat irregular since 
several of these systems are composed by a few filaments or pan- 
cakes of irregular shapes, joined by one or two vertexes. Neverthe- 
less, in order to provide a rough characterization of FVS shapes, 
we use the usually adopted ellipsoidal model to study the distri- 
bution of axes ratios. We have computed the moment of inertia 
h] = Zg.vi(^' - Cf''*') * (r^' - Cy^), where xf ' is the i-th co- 
ordinate (i.e., the x, y, or z axis) of the member galaxies, and 
are the coordinates of the geometrical centre of each FVS. The di- 
agonalization of this matrix gives the values of the three axis of the 
ellipsoid with the same moment of inertia , a, b and c. In the in- 
set of Fig. |6] we show the results of the axial ratios of FVS in an 
clb vs hi a diagram. As can be seen in this figure, neither spheri- 
cal (a ~ h ~ c) nor planar (b ~ a > c) structures are common; 
there is a preference for prolate shaped FVS structures. We use 
the triaxialit y parameter T = (I - (b/a)-)/(l - (c/df ) adopted by 
IWrav et ^ j2006i) , who analysed supercluster shapes in a ACDM 
model simulation. This parameter takes values between and 1; 
r ~ indicates oblate structures, while T ~ I indicates pro- 



late structures. Intermediate values of T represent triaxial systems, 
with triaxiality increasing (not linearly) with T. We show in Fig. 
[6] the distribution of T for samples S2, Mzs,, and M/jj,,. As it can 
be seen from this figure, the distribution of shapes is not strongly 
affected by peculiar velocities. Also, it can be appreciated that the 
observed SDSS FVS shapes are consistent with those of ACDM 
structures. The mean value of this parameter is T = 0.75 + 0.23 
for the main sample of FVS of the SDSS, and T = 0.77 ± 0.23 for 
the sample of the real-space simulated catalogue. This indicates 
that these structures are mainly triaxial systems, and the distribu- 
tion of T shows the predominance of prolate systems over oblate 
systems, in agreement with Einasto et al. (2007). The mean values 
of T reported bv lWravetalJ l i2006h range between T = 0.65 and 
T = 0. 69 depending on the linkin g length, in agreement with our 
results. ICosta-Duarte et al] ( 1201 Ol) present evidence of a trend be- 
tween supercluster luminosities and shapes, where filaments are on 
average more luminous than pancakes. We also explored the depen- 
dence of the T parameter on FVS luminosity and volume, finding a 
similar tendency. Larger values of T correspond to more luminous, 
larger, FVS. 

4.2 Testing the reliability of the method 

4.2.1 Redshift dependence 

We have also searched for a possible redshift dependence of our 
results; notice that given the shallow depth of our samples this con- 
stitutes a test of our method rather than a search for true time evolu- 
tion. To this aim, we use the samples SI and S3, defined in section 
[2l The similarity of the three histograms in Fig. |5(a)| and Fig. |5(b)| 
corresponding to these diff'erent subsamples shows a lack of red- 
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Figure 6. Probability density distribution estimates of the shape indicator 
parameter T (defined in section Section |4J). The distributions con'espond 
to FVS from the SDSS galaxy catalogue (solid line), the real-space mock 
catalogue (short dashed line) and the redshift-space mock catalogue (long 
dashed line). The inset shows a scatter-plot of the semi-axis ratios cjb and 
bla that characterize the shapes of the FVS. The sizes of circles are propor- 
tional to the number of galaxies contained in each FVS. 



shift dependence. This is also seen in the insets of Fig. |5(a)| and 
Fig. |5(b)| in the dependence of the median luminosity and volume 
with redshift. As it can be seen, there is no significant dependence 
of the median luminosities and volumes with redshift, indicating 
the reliability of the method and its lack of a redshift bias. 



4.2.2 Luminosity cut— ojf dependence 

With the aim of analyzing the effects of using tracers of differ- 
ent luminosities, we applied three magnitude cuts, corresponding 
to those of samples S2 and S3, to galaxies in the redshift range 
0.04 < z < 0.10. This defines two new samples: S2c and S3c, as 
described in Table[T] In figures [7(a)] and [7(b)] we show the distribu- 
tion of luminosities and volumes of FVS selected from the Si, S2c 
and S3c samples. We recall that we have applied the correspond- 
ing completeness luminosity factor F (Eq.[6]l to the luminosities of 
FVS. As can be seen in these figures, there is a very good agree- 
ment between the FVS derived from these different galaxy tracers. 
These tests give confidence that the use of volume limited samples 
of high luminosity galaxies do not change significantly the results, 
with the advantage that they can trace larger volumes. 



5 CORRESPONDENCE OF KNOWN SUPERCLUSTERS 
AND FVS 

As mentioned in Section [3] superclusters have been previously 
identified either as luminosity density enhancements, or as systems 
of galaxy clusters. None of the criteria ensure that the structures 
will evolve into gravitationally virialized structures. 



Therefore, it cannot be ensured that all known superclus- 
ters will be 'island universes' in the future. Rece nt superclus- 
ter catalogues were co mpiled from the 2dFGRS l lEinasto et alj 
2007). the SDSS-DR4 jSinasto et al.i ,200 6) and the SDSS-DR7 
dCosta-Duarte et al]|201(]|) . Since some of these catalogues cover 
only partially the volume of the central zo ne of SDSS-DR7, w e 
also used the supercluster identification of lEinasto et al.l ( 1200 lb . 
where the identification is based on X-ray selected and Abell clus- 
ters. Combining data from these catalogues we can correlate our 
results with previous identifications and explore which known su- 
perclusters will also be future virialized structures, at least within 
the ACDM scenario. In Fig. |5]we show the FVS associated to the 
Sloan Great Wall, Ursa Major, Corona Borealis and Bootes super- 
clusters. This figure shows the galaxies belonging to FVS asso- 
ciated to each of the above mentioned superclusters, and also the 
galaxies in neighboring FVS. For comparison, the locations of the 
centres of superclusters in the same region obtained from the liter- 
ature are presented. 

Sloan Great Wall. There are several id entified structures in 
the region associated to the Sloan Great Wall. lEinasto et all ( 1200 ih 
claim that the structures in their catalogue named SCLlll (a.k.a. 
Virgo-Coma) and SCL126 are part of the Sloan Great Wall. Of 
these superclusters, SCL126 is the richest in the Sloan Great Wall. 
lEinasto et alj (I2OIOI) find that th e structure previously identified as 
SCL136in lEinastoetal.l ( 1200 ih is in fact a part of the SCL 1 26 ,and 
SCLlll comprises three concentrations of rich clusters connected 
by filaments of galaxies. . Eina sto et aL ( 2010 ) also find the presence 
of substructures, indicating possible mergers and infall in different 
superclusters belonging to the Sloan Great Wall region, suggesting 
that the cores of these superclusters are not virialized and still as- 
sembling. We find that, according to the catalogue of FVS, SCL 111 
and SCL 126 may evolve into a single virialized structure, while the 
neighboring SCL91 could remain isolated. The structure SCLIOO 
lEinasto et alj ( l200ll) . associated with the Leo A supercluster and 
located near the Sloan Great Wall, has no FVS that it can be di- 
rectly associated with. We find several FVS surrounding the centre 
of SCLIOO thus suggesting that this supercluster could be disrupted 
in the future. 

Ursa Major Supercluster. Ursa Major is a nearby z - 0.06 
and relatively isolated supercluster. It has been found to be com- 
posed of th ree large filaments with mea n redshifts z=0 . 051, .060 
and 0.071 dKopvlova & Kopvlovl I2OO6I) . lEinasto et al] ( 1200 ll) as- 
sociate Ursa Major to the SLC109 superstructure in their cata- 
logue, and present its geometrical centre located at tt = 177.1 deg, 
S = +55.0 deg and z = 0.06. Einasto et a l. ( 2006) identifies the 
more distant filament as another individual supercluster (SCL384). 
We found four different FVS with filamentary structure associated 
to this supercluster. This indicates that the filaments would eventu- 
ally evolve into separate structures. 

Corona Bore alis Supercluster. T he Corona Borealis super- 
cluster (SCL158 in lEinasto et alj ( l200lh ) is a prominent example of 
a supercluster in the northern sky. This structure comprises 500 
galaxies at z = 0.07, with a nearly spherical morphology. We find 
that Corona Borealis is part of a larger FVS, which also includes 
the superclusters SCL805 and SCL761 identified bv lEinasto et al.l 
(2006). According to our analysis, since the total luminosity of this 
FVS is well beyond the threshold luminosity calibrated using the 
simulation, these structures are candidates to merge and form a sin- 
gle virialized system in the future. 

Bootes Superc lusters. The Bootes supercluster (SCL 138 in 
lEinasto et al] ( 1200 ll) ) is located at z ^ 0.065, and Bootes-A 
(SCL150) lies directly behind at z ^ 0.1 1. These superclusters have 
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Figure 7. Luminosity and volume distributions of structures for samples SI (dashed line), S2c (solid line) and S3c (dotted line). These three samples are 
liinited to z < 0. 1, but comprise galaxies with different absolute Hmiting magnitudes (see Table[TJ. Error bars are Poisson uncertainties within each bin. 



not been widely studied. We find that a single FVS of filamentary 
shape can be associated to each one of these superclusters. 

The comparison presented in this section helps visualize the 
diff'erences that can be found between different superstructure cat- 
alogues; our approach has the relative advantage of allowing an 
interpretation of known structures in terms of whether they will 
evolve in the future into single virialized structures. However, we 
remind the reader that this interpretation depends on the assumed 
cosmology, which in the present case is the concordance ACDM 
model as well as the several hypothesis adopted. 



6 CLUSTERS AND GROUPS WITHIN FVS 

Although it has been claimed that the clustering properties of dark 
matter haloes depend exclusively on their total mass, recent analy- 
ses of numerical simulations have shown that they can be strongly 
affected by their assembly history (see e.g. lLacerna & Padill j20ia. 
and references therein). If galaxy groups and clusters within su- 
perstructures have undergone a different evolution, it would be ex- 
pected systematic differences are to be expected at in the present 
time. The following analysis could help to address the role of the 
large-scale structure in the formation and evolution of galaxies. 
In this section we derive a simple statistical analysis concerning 
groups of galaxies and the FVS. To this end, we use three different 
samples of galaxy systems: (i) Groups of galaxies identified in the 
SDSS-DR7 (see Section[2l(, suitable to search for correlations be- 
tween halo mass and the large scale environment they inhabit, (ii) 
Abell AGO cluster catalog ( Andernach 1991), that provides redshift 
measurements of 1059 clusters in the footprint area of the SDSS- 
DRV, and (iii) X-ray clusters of galaxies, from the RASS catalogue 
jPopesso et alj|2004h providing a suitable sample of bright X-ray 
clusters within the SDSS area. 



We have computed the fraction of galaxy groups/clusters 
within FVS as a function of their mass (Fig.|9](. In this figure, error 
bars indicate Poisson uncertainties given by, 

^(M^^r) = V<^/^"' (7) 

where A',',','^ is the number of groups/clusters inside a FVS in the 
i-th mass bin, and A'*',] is the total number of groups in the same 
mass bin. As can be appreciated, the fraction of groups within FVS 
is greater for higher group masses. This result is also present in the 
catalogue of Abell-AGO clusters. This catalogue supplies 177 clus- 
ters within the coverage area of SDSS-D R7, out of which 11 9 lie 
within the supercluster sample defined by Einasto et al.l ( 1200 ih and 
63 lie within FVS. A total number of 55 Abell clusters are found in 
both catalogues. We compute fraction of Abell clusters within FVS 
as a function of the richness class parameter. We find only 14.5 per 
cent of type clusters within FVS; a value that increases for richer 
Abell clusters: 24.4 per cent for type 1, 34.4 per cent for type 2 and 
41.7 per cent for type 3 clusters. Thus, there is a significantly larger 
probability of the richest Abell clusters to be in FVS, in contrast to 
poor clusters. 

It is also worth studying X-ray emitting clusters of galaxies. 
While Abell clusters have an estimate of their richness that is dis- 
crete and approximate, the X-ray cluster catalogue provides X-ray 
luminosities, known t o be good indicator s of the underlying gravi- 
tational potential well jRvkoff et al.l200^ . The catalogue of X-ray 
clusters also allows to study the fraction of clusters inside FVS as 
a function of the X-ray luminosity. We obtain the same trend in 
this fraction as in the case of DR7 groups, as it can be seen in Fig. 
[Tol However, while the fraction of systems inside FVS is at most 
0.3 for the DR7 groups (except for the last virial mass bin, which 
has a large uncertainty), we find that 85 per cent of the clusters 
with Lx > 4 10^' W are part of a FVS. A similar value is found for 
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Figure 8. Spatial distributions of galaxies within Future Virialized Systems associated to known superclusters (filled circles). The positions are shown in a 
cartesian system associated to the equatorial coordinate system (x-axis pointing to the vernal point y and z-a xis pointing to north equatorial pole); the line 
indica tes the observer direction. We also show the centres of superclusters in the catalogues of Ein asto et alj ^2001). Einasto et al. (2007) and Einas to et alj 
Galaxies within nearby FVS are indicated with black dots. In the Sloan Great Wall region (a), we identify a single FVS associated with SCLlll, 
SCL126 a nd SCL136, indicating that t hese superclusters may evolve into a single viriaHzed system. The Ursa Major Supercluster, in (b), consists in three large 
filaments JKodvI ova & KoDvlovl200^ . However, we find a different FVS associated to each filament, suggesting that the filaments would evolve into separate 
structures. In (c), the Corona Borealis Supercluster SCL158, is part of a larger FVS, that also includes SCL805 and SCL761. So, these three systems could 
merge to form a single virialized structure in the future. Bootes and BootesA Superclusters in (d) have one-to-one correspondence with FVS, so each system 
would evolve into a future virialized structure. 
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Figure 9. Percentage of DR7 galaxy groups inside FVS witli respect to the 
total number of groups into virial mass bins. The first and second mass bins 
comprise 21 galaxy groups, although none of them are located into a Future 
Virialized Structure. 

superclusters by lEinasto et alj ( 1200 ih using a compilation of X- 
ray clusters based on the ROSAT All Sky Survey. The authors also 
found that in rich and very rich superclusters the fraction of X-ray 
clusters is higher than the fraction of Abell clusters, consistent with 
our findings. 

The results of this section are in agreement with the expecta- 
tion of an assembly bias scenario, where the oldest, most massive 
systems are preferentially located in present-day overdense regions. 
Also, this suggests that the most frequent future cluster merger 
events will be associated to the most massive clusters within FVS. 



7 DISCUSSION 

On the observational side, catalogues of superstructures have 
been presented on increasingly larger galaxy datasets (see 
iLiivamagi et all l2010l and references therein). These catalogues 
provide important characterizations of the largest scale structures 
and allow for different tests of structure formation models. Regard- 
ing the evolution of these structures from the present, numerical 
simulations reveal the conditions by which these systems would 
evolve in to isolation and dyn amical equilibrium within the ACDM 
scenario. iBusha et alj j2005l) consider the long term evolution of 
large structures, up to a = 100, using numerical simulations, and 
analyse the definition of a radius that encloses all the mass that ulti- 
mately will form an isolated structure. The authors find that a tran- 
sition region is found during the matter dominated era (a < a^,), 
between an inner hydrostatic region and an outer region that ex- 
pands with the perturbed Hubble flow. In this intermediate region, 
accretion of matter towards the mass concentration is complex and 
gives rise to a variety of definitions of scales to describe clusters. 
At later times, the accretion region narrows and the hydrostatic 
and turnaround scales merge to form a single zero-velocity sur- 



Figure 10. Percentage of X-ray clusters in FVS with respect to the total 
number of X-ray clusters, as a function of X-ray luminosity. 

face th at unambiguously defines the halo mass. lArava-Melo et al.l 
( l2009l) argue that, once A has started to dominate the expansion of 
the universe, the cosmic web growth stops, so that the spatial dis- 
tribution of superclusters is essentially the same starting from the 
present epoch. The authors identify sup erclusters on the basis of 
criteria presented in lPiinner et al ] ( l2006h . These criteria, based on 
the spherical collapse model and implemented on N-body simula- 
tions, allow to develop a method to isolate structures in the present 
universe in redshift-space that will eventually f orm a bound struc- 
ture. In a subsequent work. lPUnner et aL presents a detailed 
method to be applied in observational data. This method is based 
on the assumption of spherical collapse and makes use of projected 
velocity envelope templates and mass profiles. 

The "gravitationally bound" cr iterion is widely used, and was 
also applied to our nearby Universe. iNagamine & Loebl ( l2003h per- 
formed a numerical simulation with initial conditions at z = re- 
producing the observed galaxy distribution in the local Universe. 
They conclude that the mass overdensity of the Local Group is 
above the required threshold to be a bound structure in the future. 
According to the authors, Andromeda and the Milky Way will prob- 
ably merge, while our Local Group and the Virgo Cluster will not 
form a virialized structure in the future^ 

From a theoretical point of vie w.lEinasto et al] ( l2010h decom- 
pose the luminosity density field of lLiivamagi et al.l ( 1201 Oh using a 
wavelet analysis. They also study the formation of the cosmic web 
based on the evolution of the density perturbation phases in numeri- 
cal simulations. They state that very rich superclusters are the result 
of large scale {X > 32h^^ Mpc) density waves, combined in similar 
phases. Also, the cosmic web develops in early stages of the Uni- 
verse evolution, and it is generated by the phase synchronization. 

The largest supercluster catalogue constructed so far with 
SDSS-DR7 data dUiivamagi et al.l l2O10h allows to test our pro- 
cedure and compare the results. As they apply the density field 
method using a wide range of threshold density parameters (ana- 
logues to our Dc), we can directly contrast the final results of both 
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Figure 11. Upper panel: luminosity (a ) and volume (b) for FVS obtained 
from the S2 sample (s olid line) and for Superclusters in the catalogue of 
iLiivamagi et al.l )2O10l) (dotted line). Error bars are Poisson uncertainties 
within each bin. Lower panel: Fraction of galaxies in Superclusters that 
also belong to a FVS, fiN. 80 per cent of superclusters have at least a 65% 
of their galaxies within FVS, represented by the vertical dashed line. 

catalo gues for the same value of = 5.5. The ILiivamagi et aEI 
( l2010h catalogue is deeper, so we restrict the analysis to z = 0.12, 
the limit of our main sample of FVS. Up to this di stance, we iden- 
tify 1 50 FVS compared to 142 superclusters in ILiivamagi et al.l 
1 I2OIO ) catalogue. In the upper panel of Fig. [TT] we show lumi- 
nosity and volume of these structures. From the available list of 
galaxies in superclusters of ILiivamagi et a i] ( |2010|) . we selected a 
subsample of 83.798 galaxies, with the same limits of our sample 
S2 (0.04 < z < 0.12 and < -20.47). We find that 79.4% of 
the galaxies in superclusters also belong to a FVS, and only the 
3.4% of galaxies outside of superclusters are associated to a FVS. 
We also analyse the galaxies within each supercluster separately, 
defining the fraction of galaxies in the supercluster that are part of 
a FVS fiN, and the fraction of galaxies in the supercluster that do 
not belong to any FVS, four- 

In the lower panel of Fig. [TT]we show the //a, fraction for 
the galaxies within each supercluster. We find that 80% of super- 
clusters have //jv over 0.65, and four lower than 0.35. The pre- 
vi ous analyses ind i cate t hat the general features of the structures 
of ILiivamagi et al.l bOld) and those of the FVS show a very good 
agreement. 



8 CONCLUSIONS 

The properties of the largest structures in the Universe provide im- 
portant information on the cosmological model. In this work we 
have developed a physically motivated procedure to isolate struc- 
tures in the nearby universe that will evolve into gravitationally 
virialized systems in the future, within a AC DM framework. The 
method starts from a smoothed luminosity density map, where 



high density peaks are separated as superstructure candidates. This 
method depends on a number of free parameters that we deter- 
mined on a theoretical basis. We adopted an Epanechnikov Ker- 
nel following recent similar implementations of the method (e.g. 
lEinasto et allbOOTl : ICosta-Duarte et al.ll2O10h . We then calibrated 
the luminosity overd ensity parameter Dj according to results from 
IPtinner et alj ( 1200^ . who used numerical simulations to establish 
the conditions for a superstructure to be virialized at a ~ 100. 
We find that a cut in the total luminosity of superstructures at 
L = lO'^Lo is convenient to additionally clean the sample of FVS 
(see Fig.|2]l- This is due to the fact that although some superstruc- 
tures have a total luminosity overdensity above the critical value, 
their mass can still be not enough to ensure the future collapse onto 
a virialized system. We used a mock catalogue to further constrain 
and test the reliability of the identification method. We defined 
three samples of FVS from volume limited samples of galaxies. 
From the comparison of these samples we conclude that the method 
is not affected by the constraints in the survey redshift and the lu- 
minosity of galaxies, since the small differences can be explained 
by the corrections implemented assuming a universal luminosity 
function. 

The joint analysis of group catalogues and the samples of FVS 
indicate that the fraction of groups within FVS increases with the 
virial mass (see Fig.|9]l- The fraction of X-ray clusters belonging 
to FVS is significantly larger, and also increases with the X-ray 
luminosity (Fig.llOt. This is consistent with the assembly bias sce- 
nario, with massive, old assembled systems located preferentially 
in present-day superstructures as FVS. Thus, future cluster merger 
events will often involve the most massive, bright X-ray clusters. 
We have also analysed particular superclusters previously studied 
by other authors. Many known superclusters are found to be coin- 
cident with one or more FVS in our sample. As an example, our 
analysis shows that the Ursa Major supercluster is coinposed by 
four filaments, each one a distinct FVS; and the Corona Borealis 
and Bootes superclusters are clearly associated with FVS in our 
catalogue. 
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